Discrete instability in nonlinear lattices 
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The discrete multiscale analysis for boundary value problems in nonlinear discrete systems leads 
to a first order, strictly discrete, modulational instability (disappearing in the continuous envelope 
limit) above a threshold amplitude for wave numbers beyond the zero of group velocity dispersion. 
Applied to the electrical lattice [Phys. Rev. E, 51, 6127 (1995)], this acurately explains the experi- 
mental instability at wave numbers beyond 1.25 rad.cell -1 . The theory is also briefly discussed for 
the sine-Gordon and Toda lattices. 
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Introduction. The Benjamin- Feir instability describes 
the mechanism of the exponential growth of the modu- 
lation of a wave evolving in a continuous nonlinear dis- 
persive medium It has revealed as a fundamen- 
tal process occuring in various physical situations, from 
water waves to nonlinear optics through plasma waves, 
electrical transmission lines, etc., see e.g. |]. This in- 
stability, originally established for periodic wave trains 
(Stokes waves in surface water theory) Q, is actually a 
generic process in physics, well understood in the context 
of the nonlinear Schrodinger equation (NLS) |5j ||. In 
the case on an infinite medium, and for the focusing NLS, 
this instability has been shown to evolve eventually into 
localized coherent nonlinear structures (the solitons), see 
e-g- §• 

The analysis of a nonlinear system is performed by first 
deriving from the original model a simpler limit equation 
which usually results to be the NLS. This universal limit 
is obtained by a very general approach called multiscale 
perturbation analysis M. It is worth remarking already 
that in such cases the instability is actually a second or- 
der effect: NLS is obtained as a first order perturbation 
and the Benjamin- Feir instability results from a pertur- 
bative analysis of NLS itself. 

In the case of discrete systems where intrinsic nonlin- 
ear modes (called also discrete breathers) are known to 
exist ||- [Q, the modulational instability, mechanism 
for their spontaneous generation out of an initial distur- 
bance, is usualy described in a semi-discrete multiscale 
analysis where the envelope of the discrete carrier wave is 
treated as a continuous function Jll|- p7|l . Fully discrete 
treatments have also been proposed as in |L8 19 1, based 
on the rotating wave approximation which essentially con- 
sists in neglecting harmonics. In a nonlinear system, this 
procedure evidently may give wrong predictions for large 
times. 

In some of these studies however, the problem is not 
that of the long time evolution of an initial perturba- 
tion but rather the far behavior of an input boundary 
disturbance (or forcing). This is particularly evident in 
the electrical transmission line of |lj|[l8| where an input 
signal (modulated wave) is sent at one end of the lattice 
and collected at various locations along the line. 

Hence we consider here the problem of the scattering of 
a plane wave incoming into a nonlinear discrete medium. 



Based on the method founded in (gCfl , a discrete change of 
variables is introduced and used to construct the related 
discrete limit model. This results as a NLS-like equation 
in discrete space and with space and time exchanged. For 
wave numbers beyond the zero of the group velocity dis- 
persion, it is shown that the waves are uniformly unsta- 
ble as soon as their amplitudes exceed a threshold value. 
Moreover this instability is a purely discrete property as 
it is shown to disappear if the envelope is treated as a 
continuous function like in the semi-discrete approach. 

This approach is applied to the electrical lattice of JlJ] 
where the problem is indeed that of the far behavior of 
an input boundary datum and for which a set of very de- 
tailed and acurate experimental results are displayed. In 
particular the measured threshold of modulational insta- 
bility occurs at wave numbers beyond 1.25 (in rad.cell -1 ) 
which is the value predicted by the present theory. We 
predict also a threshold amplitude depending on the in- 
put carrier frequency and modulation. For one of the 
experiments of this threshold value is 0.23 V, signif- 
icantly smaller than the amplitudes actually used in the 
experiments where damping had to be overcome. This 
threshold amplitude is shown to depend on the relative 
frequencies of the carrier v.s. modulation and hence can 
be easily varied which should motivate new experiments. 

Basic tools. Treating a boundary value problem cor- 
responds to making an expansion of the wave number 
around small deviations of the frequency from the linear 
dispersion law. To see this, we consider a plane wave 



A(n,t) = exp[i(Kn - fit)] 



(1) 



in a nonlinear lattice (with intcrsitc spacing normalized 
to 1) with the linear dispersion relation of bandpass filter 
type 



ti 2 



wl + 4?^ sin 2 (if/2) . 



(2) 



Due to nonlincarity, the actual frequency to and wave 
number k experience deviations form the values f2 and 
K obeying the linear dispersion relation (||). There are 
two ways of representing these deviations. 
With the wave packet 



,(*) 



dk u{k)e i( - kn - ut \ 



(3) 
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we may expand the frequency in terms of the assumed 
small deviations from the value K of the wave number as 



k = K + e\ , ui = SI + ev„ X + e 2 P A 2 



(4) 



where v g — dVl/dK is the group velocity and 2P = 
d 2 il/dK 2 the group velocity dispersion. In that case, 
the expression for the wave packet becomes 

u n (t) = ee^ Kn - nt U(Z,T), (5) 

where the slowly varying envelope <p, namely 

<j>(Z,T) = [ d\u(X)e^ xz - x2pT ^ (6) 



(7) 



is written in the frame 



Z = e(n — v g t) , T = e 2 t. 



This means that the effects of an initial disturbance are 
observed for large (e -2 ) times by following the waves at 
the group velocity and analyzing the perturbation about 
the plane wave solution. This approach then corresponds 
to treating initial value problems. Note in particular that 
the variable Z hereabove defined becomes by construc- 
tion a continuous variable. 

Alternatively, considering the wave packet 



u n {t) = / du uHe 1 ^"-^, 



(8) 



we may expand here the wave number in terms of the 
frequency as 



lu = £1 + ev . k = K + e h e 2 Q v 1 

Vn 



(9) 



(with 2Q = d 2 K/dfl 2 ), hence giving 

u n {t) = ee* Kn - m H{tr), (10) 
with the slow modulation 

ip(Z,r) = [ du u{v)e i{ - v * Q t- VT \ (11) 



in the frame 



e(i - n/vg) , i = e 2 n. 



(12) 



Such a reference frame has a clear physical meaning: one 
considers long distance (e -2 ) effects in the retarded time 
to let the input disturbance enough time to reach the 
observed lattice point. Hence this corresponds to the sit- 
uation where the lattice is excited at one end, in other 
words to a boundary value problem. This is the tool we 
shall use hereafter when the variable £ is kept discrete as 
follows. 



Discrete multiscaling. Performing a discrete multi- 
scale analysis goes first with assuming a hypothesis of 
slow modulation. This means for the envelope ip n (t) of 
the wave A(n,t), that 



||9^„(i)||~e||cU(M)|| 



(13) 



which then allows to chracterize the dimension e of the 
slow variable r by ( p"2| ) . 

Second, as introduced in |^p|] , we define a large grid in- 
dexed with the slow variable m defined by sampling the 
original grid at each N = e~ 2 (more precisely the closest 
odd integer) . As a consequence we can index the variable 
£„ by m in the new grid, namely we can call m a given 
point n and m + j the points n + jN for all j. In short 



1, £n = m, t; n+N = m + 1, 



(14) 



Considering now a slow modulation ip(^ n ,T n ) in the 
variables (12), of the plane wave A(n,t), we represent 
this modulation by the function if>{m, r) of the above 
defined discrete variable m and of r = t„. Then it is 
demonstrated in ||^] that the first centered difference of 
u n {t) = A(n,t)i/j(m,T) reads 

Un+l-Un-l = [A n +1 - A n -i}lp m 

+ {e/v g )[A n+1 + A n -i]d T ip m 



+e 2 /2[A n+ i + A n ^x]bPm+l ~ tpm-l] 



(15) 



while the second discrete derivative is given by 



u n+l -2u n + u„_x = [A n+1 - 2A n + A n _x] -0m 

+ {t/Vg)[A n+1 - A^^drlpm 

+(t/v g ) 2 ^[A n+1 + A„_i]d 2 m 
+e 2 /2[A„ +1 -A„_ 1 ][^ m+1 -Vm-i]+0(e 3 ) • (16) 

The discrete electrical lattice. We consider now the 
problem of the onset of modulational instability in the 
nonlinear electrical lattice experimentaly and theoreti- 
caly studied in iQ. The model reads 



(A + V n )d 2 V n - {d t V n f 



A 



v, 



n+1 



+ V n -! - (2 + 



(17) 



where A is a constant and V n (t) is the voltage of the n-th 
cell. The linearized version of this equation does possess 
the dispersion relation (^). The experiments consist in 
sending at one end of the electrical lattice a voltage V (t) 
as a slightly modulated plane wave for various values of 
the carrier frequency and amplitude, and then measuring 
the output voltage at some cells along the lattice. 

To proceed with discrete multiscaling, it is more con- 
venient to work with the equivalent system 
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d t V n =u [l+ \V n ] B n , 



d t B n = u 



K.+i + K-i - (2 + %)V n 



(18) 



A perturbation expansion is now performed with the fol- 
lowing series 



(19) 



where of course A (n,t) — exp[i£(Kn — fii)], and where 
the reality condition results in ipp ^ = ipp and <j)p ^ — 

Inserting the above expressions in ( p"8| ) and using the 
formula (16), we arrive at 
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(20) 



where r\ = rj m (T) and 5 = 5 m {r) are the new unknowns 
and where we have defined 

a = -3cjq - 4uq(1 - cos if) 2 . 

Solving then the system for <f>^ and t/^ , the unknown 
function <5 m (r) cancels out and we are left with the equa- 
tion for rj m (r) which eventually writes 



-[rjm+l - Vm-l] + Q -q^T - 7\Vm\ Vm = 



with the following if -dependent coefficients 
_ 1 d 2 K _ itg(l - cos Kf - uj 2 cos K 



(21) 



2 dfl 2 
1 



2u* sin 3 K 
fl 2 [u 2 + 4u 2 sin 2 K] 



' 2A 2 u 2 smK[3Lu 2 + 4u 2 (l-cosK) 2 } 
From ( |l9| ) and (|2(]) we have then 

V n (t) = e v m (r)e^ Kn ^ + 0(e 2 ) + c.c. 



(22) 



(23) 



Note that higher order corrections to V(n, t) are either 



explicitely expressed in terms of rj through (20) or by 
linear inhomogeneous differential equations. That means 
that the theory is self-consistent and in particular that 
the overtones are by no means neglected. 



Modulational instability. The continuous version of 
the equation ( pl| ) is a well known model for boundary 
value problems in optical fibers |2l|j . It has been stud- 
ied also in the case of initial value problems in unstable 
media to describe for instance Rayleigh- Taylor insta- 
bility. 

The plane wave solution to (|2l|), say 



?7 m (r) = Bexp[z(Am - ut)] , 
obeys the dispersion law 

i , 2 = -I[sinA + 7 |i?| 2 ] • 



(24) 



(25) 



As we are considering a boundary value problem, v and 
B are the given frequency and amplitude of the modu- 
lation. Then, as sin(A) is bounded, and this is precisely 
a discretreness effect (see conclusion), there is no real 
solution A if 



Q > 



\B\ 2 > - = B 2 cr , 

7 



(26) 



leading to unstable plane waves. 

Reporting to the definition ( p2] ) we first remark that Q 
vanishes at K = K cr which is also the zero of the group 
velocity dispersion. Second, it is easy to see that 



K < K cr =^ Q < , K > K cr Q > 



(27) 



Then the waves are unstable for K > K cr and \B\ > 
\B cr \. 

Application. The parameters of the experiments of 
@ are 

lu = 2.13 MHz , u = 1.77 MHz , A = 3.9 V , (28) 
for which the zero K cr of a(K) results to be 

K cr = 1-25 rad.ceir 1 , (29) 

in agreement with the experimental observation of |l8| . 

Concerning the threshold amplitude B cr defined in 
(^6|), its actual value in the physical units is found by 
going from rj m (r) to the value of V n (t) through (pq), 
namely the critical value of the amplitude of V n (t) is 
simply eB cr . Hence we must evaluate e by using its very 
definition, the requirement of slow modulation (|l3|). In 
one of the the experiments of (see their fig. 6), an in- 
put plane wave of frequency / = 600 Khz (wave number 
K e = 2.15 rad.cell -1 ) is modulated at F = 32 Khz and 
sent in the lattice. Remembering the change of variables 
(|l2|), in particular r = et, such a modulation implies 
e = f/F ~ 0.053. 

For the parameters (|28]), the critical amplitude B cr , at 
wave number K e is easily evaluated and hence value of 
the instability threshold results to be 



eB cr (K e ) = 0.23 V 



(30) 



This value is significantly smaller than those currently 
used in the experiments of uM where damping prevents 
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the use of smaller voltage input. However, the actual de- 
pendence of B cr on the relative frequencies of the carrier 
and modulation waves (through the parameter e) allows 
one to vary the value of threshold B cr . This property 
could be used to check the threshold amplitude on new 
experiments. As shown on the figure, we note finally that 
the function B cr (K) is bounded (maximum of 4.68 V at 
K = 2.59), hence this threshold of instability exists at 
any carrier wave number in ]K cr , tt[. 
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Fig.l: Critical amplitude B cr (K) in units of e 
Predictions. The theory can also be applied to the 
sine- Gordon chain 

u n - Mq (u n+1 - 2u n + u n -i) + sinii„ = , (31) 

which has the linear dispersion relation (||) and for which 
the limit model is (21) with the coefficient 



7 : 



4uq sin K 



(32) 



and of course the same value for Q. Performing the anal- 
ysis shows an instability for wave numbers greater than 
the zero K cr of Q(K), which develops for modulation am- 
plitude obeying (in units of e and in the Brillouin zone 
Ke}0,n[) " 



B > B rr = 



2uq^so\.K 



For the integrable Toda lattice whose dispersion rela- 
tion is 17 = 2sin(if/2), leading to the limit model ^T" 
with (20 1 



Q 



1 sm{K/2) 
' 4 cos 3 (K /2) ' 



1 



7 : 



sin if 



(34) 



the theory implies that the waves are unstable for all K 
and \B\ 2 > \wa.K\. 

Conclusion. The instability described above results 
from wave scattering in a nonlinear medium (boundary 
value problem) which requires to use a change of vari- 
ables like (Of) which in turn leads to the unstable wave 



equation ( pl| ) for the envelope. The resulting instability 
is then of first order (a second perturbation analysis of 
the equation for the envelope is not required), and has 
for boundary value problems the same universal charac- 
ter as the Benjamin-Feir instability has for initial value 
problems. 

Moreover this instability is a purely discrete phe- 
nomenon as indeed, performing a semi-discrete (discrete 
carrier and continuous envelope) analysis with the new 
variables (|l^) leads to the continous version of ( pl| ) and 
hence to the dispersion law 



= > 



(35) 



The point is that, even for positive values of Q, there 
always exists a real solution A of the above dispersion 
relation for any given modulation frequency v. 
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